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Abstract 

We consider a dissipative vector field which is represented by a nearly-integrable 
Hamiltonian flow to which a non symplectic force is added, so that the phase space vol- 
ume is not preserved. The vector field depends upon two parameters, namely the per- 
turbing and dissipative parameters, and by a drift function. We study the general case 
of an ^-dimensional, time-dependent vector field. Assuming to start with non-resonant 
initial conditions, we prove the stability of the variables which are actions of the conserva- 
tive system (namely, when the dissipative parameter is set to zero) for exponentially long 
times. In order to construct the normal form, a suitable choice of the drift function must 
be performed. We also provide some simple examples in which we construct explicitly 
the normal form, we make a comparison with a numerical integration and we compute 
theoretical bounds on the parameters as well as we give explicit stability estimates. 
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1 Introduction 

A breakthrough in the theory of the stability of nearly-integrable Hamiltonian systems was 
achieved by the seminal works of A.N. Kolmogorov ([32]), V.I. Arnold ([1], [2]), J. Moser 
([35]) and N.N. Nekhoroshev ([36], [37]). Nearly-integrable systems can be modeled by 
Hamiltonian functions of the form 

H{y,x,t)^Ho{y) + eHi{y,x,t), (1) 

where, for an ^-dimensional system, y e M^, (x, t) G T^+\ Hq, Hi are regular functions. Hi 
being periodic in x and t, and £ is a small parameter, called the perturbing parameter. As far as 
£ = 0, Hamilton's equations associated to (1) are integrable. The actions are constants and the 
motion is described by periodic or quasi-periodic solutions, according to the rational or irra- 
tional character of the frequency vector, say u}{y) = dHo{y) / dy. In this context, under general 
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assumptions, KAM theory yields the persistence of invariant tori for sufficiently small values 
of the perturbing parameter, which implies a stability property for low-dimensional systems, 
while higher dimensional models may admit diffusion through the invariant tori. A stability 
result for arbitrary dimensions can be obtained through Nekhoroshev's theorem, which guar- 
antees the stability of the actions (namely the confinement in a given domain of the action 
space) for exponentially long times. Both KAM and Nekhoroshev's theorems are constructive, 
in the sense that they can be explicitly implemented to provide bounds on the parameters en- 
suring the persistence of invariant tori or the confinement of the actions over exponential times. 
For this reason the theorems have been widely used to investigate physical systems in differ- 
ent contexts. Most notably. Celestial Mechanics was a spur for the development of analytical 
results about the stability of nearly-integrable systems. In this field, there exist many appli- 
cations of KAM theorem (compare with [12] and references therein) and of Nekhoroshev's 
theorem (see, e.g., [4], [6], [7], [15], [16], [23] [27], [28], [33]); most of the latter papers deal 
with the stability of the triangular Lagrangian points. All these results are obtained in a conser- 
vative framework. However, as it is well known, many interesting physical systems of Celestial 
Mechanics are affected by a small dissipation. We quote, for example, the three-body problem 
with Poynting-Robertson drag ([3], [20]) or the spin-orbit model (compare with [34]) with 
tidal torque ([12]). These problems can be modeled by a nearly-integrable system to which a 
small dissipation is added (see, e.g., [9]). 

While in the conservative setting we find periodic orbits, invariant tori as well as chaotic mo- 
tions, in the dissipative context (where the phase space volume is not preserved by time evo- 
lution) we speak of periodic attractors, quasi-periodic attractors and strange attractors. Beside 
the case of contracting and expanding systems, we consider also oscillating systems for which 
the energy varies periodically around a mean value. Coexisting attractors have been established 
in, e.g., [25], [12]; periodic attractors are shown to exist within parameter regions known as 
Arnold's tongues ([9], see also [14]). Concerning invariant attractors, KAM results (both ana- 
lytical and numerical) provide their persistence under general assumptions (see [8], [10], [11], 
[13]); an application of the converse KAM theorem (about the non-existence of quasi-periodic 
attractors) has been presented in [19]. 

In this paper we exploit exponential type estimates for nearly-integrable dissipative vector 
fields, which are defined as follows: we consider a vector field depending on two parameters, 
say e and /x, such that for /x = we obtain a Hamiltonian vector field in action-angle variables, 
say {y, x) e X T^; we consider the most general case of a vector field depending also 
periodically on the time t. We assume that for // 7^ the vector field is not conservative; at 
the first order in /i the vector field contains a drift function, say 77 = rj{y, x, t), which must 
be suitably chosen in order to meet some compatibility conditions, which allow to construct a 
proper normal form. In this context we establish stability estimates by proving the confinement 
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of the actions of the integrable approximation within a given domain of the phase space for 
exponentially long times (see [5], [26] and especially [38] which is at the basis of the present 
work) in the non-conservative non-integrable system. The proof relies on the construction of 
a proper normal form, which results from the composition of a conservative transformation 
(removing the conservative terms to suitable orders in e) and a dissipative transformation, 
which acts on the /x-dependent terms. Perturbative methods for vector fields have also been 
developed in [24] (see also [21], [30], [31]). The case of resonant initial conditions is treated 
in [17], where one needs to construct a Lyapunov function (i.e. the energy function associated 
to the conservative system with = 0), which must be used in order to bound the variation 
of the normal form coordinates; in order to achieve the result, the resonant case requires to 
work in the extended phase space and to impose the quasi-convexity assumption, which is not 
needed in the non-resonant case treated here. The final results are conceptually different, since 
in the non-resonance case we obtain a stability result valid for exponential times, while in the 
resonant case of [17] the stability time may vary linearly or exponentially with the parameters, 
according to the structure of the vector field and to the choice of the resonance condition. 
We believe that an explicit construction of the normal form and of the stability result in the 
non-resonant case is certainly of great interest in view of applications to concrete models. 
Just to quote a field which is familiar to the authors, a non-resonant application in Celestial 
Mechanics can be of interest in order to bound the motion of the majority of the main belt 
asteroids. 

In this paper we provide also concrete examples by investigating two model problems: a strictly 
dissipative vector field and a system with oscillating energy. We test the accuracy of the nor- 
mal form construction by comparing the analytical expression with a numerical integration 
(see [18] for a discussions about the integration algorithms in nearly-Hamiltonian systems). 
Furthermore, we compute concrete estimates on the parameters (namely, the radius of the ac- 
tion domain and the stability time) and we investigate their dependence on quantities like the 
normal form order, obtaining stability times which grow exponentially with the normalization 
order. 

This paper is organized as follows. In Section 2 we set-up the definitions and notations. In 
Section 3 we state the normal form Lemma and the main Theorem; the proofs are given in 
Section 4. An application of the normal form Lemma and the computation of the stability 
estimates according to the results of the Theorem are given in Section 5 . 

2 Set-up: notations and assumptions 

We consider the following ^-dimensional, time-dependent vector field 

X = uj{y) + ehio,y{y,x,t)+ iifoi{y,x,t) 
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y = -ehio,:,{y,x,t)+iJ,[goi{y,x,t)-r]{y,x,t) 



(2) 



where y E A with A being an open domain of M^, {x,t) G T^+^, ^ g M+, /i G M\{0}, 
and 1] are real-analytic, f-dimensional vector functions with components {uj^^\ and 
{ri^^\ ...,r]^^^). We assume that hio, /oi, fifoi are known periodic, real-analytic, ^-dimensional 
functions defined on A x T^+^. In all this paper we adopt the following notations and defini- 
tions. 

(i) The subscripts x, y denote derivatives with respect to x, y. 

[ii) For integers j, m, the symbol Fjm denotes that the function Fj^ is of order f/^. 
{Hi) We say that a function F is of order km e and jji, in symbols F G Ofc(e, jj), if its 

Taylor series expansion in e, pL contains only powers of e^f/^ with j + m > k. 

(iv) For any integer vector k — {ki, ...,ke) e we introduce the norm |A;| = + ... + 

(f ) We denote by a bar the average of a function over the angle variables, while the tilde 
denotes the oscillatory part; more specifically, we decompose a function F = F(y,x,t) as 

F{y,x,t)^F{y) + F{y,x,t), 

where the average F(y) is given by 

F{y) = {F{y, X, t))^,t = TTTwr / ^' ^) ^^^^ ' 

while the oscillatory part is defined as F{y, x, t) = F{y, x, t) — F{y). 

We assume that there exists a subset D Q A such that the vector function u — ui{y) satisfies 
the following non-resonance condition up to a suitable order K with X G Z+: 

\oj{y) ■ k + m\ > a ioraWyeD, A; G , meZ, |A;| + |m| < i^T , (3) 

where a is a strictly positive real constant and the dot denotes the scalar product. 

With reference to (2) we call e the perturbing parameter, while we refer to /x as the dissipative 

parameter For e = = the equations (2) are trivially integrated as 

x{t) = x{0) + uj{y{0))t 

y{t) = y{o), 

where x{0),y{0) are the initial conditions at time t = 0. For e small and = 0, the equa- 
tions (2) reduce to the conservative vector field associated to the nearly-integrable Hamiltonian 
function 

H{y, X, t) = hoo{y) + ehio{y, x, t) , 
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where /iqo is such that uj{y) = /ioo.y(|/)- For 7^ we assume that the vector field is dissipative 
for any y E A, (x.t) E namely there exists a subset of the phase space, which is 

contracted or expanded asymptotically by time evolution into a compact set ([39]). We also 
consider systems such that the energy is oscillating around a mean value. 

We refer to 77 = r]{y, x, t) as the dWft vector function with components {r]^^\y, x, t), ri^^\y., x, t)), 
which can be expanded as 



00 m 



m=0 j=0 

We remark that is an unknown function, which will be properly chosen in order to meet some 
compatibility requirements in order to perform a suitable normal form (compare with KAM 
results like in [11]). 

Let a function / = f{y, x, t) be defined for y E A, (x, t) E T^^^; we denote by Cro{A) the 
complex neighborhood of A of radius tq, namely 

Cro{A) = {yEC^: \\y - yA\\ < ro for all yA & A} , 

where || • || denotes the Euclidean norm. Let Cgo (T^"*"^) be the complex strip of radius sq around 
T^+^, namely 

C,,{T'+') = {{x,t) E C'+' : max \^{xj)\ < Sq , |»(t)| < Sq} , 

i<j<e 

where denotes the imaginary part. Let us denote the Fourier expansion of a function / = 

f(y,x,t) as 

f{y,x,t)^ J2 A^(y)e^^'-«^ 

(fc,m)eZ^+i 

For an analytic function on A x T^+^ we introduce the norm 

sup Yl lA-(?/)|e(l^l+H)^°, 

while for a function g = g{y), g : M.^ ^ R, we define lls'llr,, = ^^Py&CrQ{A) \\9{y)\\- For an 
^-dimensional vector function v with components (vi, v^), we define 



Ell". 
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For a function / = /(y, x, t), / : x T^"*"^ — )■ M, and for any positive integer K, we denote 
by f-^, f^^ the sum over the components with Fourier modes less or equal, or respectively 
greater than K, namely 

(fe,m)eZ^+i, \k\+\m\<K 

and 

(fe,m)eZ^+i, \k\+\m\>K 

3 Stability for exponential times: statement of the results 

Stability estimates for exponential times are obtained by implementing a change of coordinates 
such that the vector field (2) is transformed to a suitable normal form of order N. Precisely, let 
us consider a coordinate transformation close to the identity, say S*^") : A x T^+^ 
such that 

(F, X, t) = s(^) {y,x,t), r e M^ {X,t)e T^+^ , (4) 

where S^^^ depends parametrically also on e, /i; we require that the transformation acts as 
the identity on the time. Let us assume that on a suitable parameter domain, the coordinate 
transformation (4) can be inverted as 

y^y{Y,X,t), x^x{Y,X,t). 

In the following Normal Form Lenama we look for a change of coordinates (4) such that (2) is 
transformed to the following normal form of order N: 

X = Qf\Y;e,fi) + FN+,{Y,X,t) + F>''iY,X,t) 

+ 5/i>^,(y(F, X, t),x{Y, X, t),t) + /x/o>i^(y(F, X, t),x{Y, X, t),t) 

Y = GN+i{Y,X,t) + G>f{Y,X,t) 

- eh>,%{y{Y, X, t),x{Y, X, t),t) + i^g>''{y{Y, X, t),x{Y, X, t),t) , (5) 

where Fjv+i, Gn+i are vector functions of order On+i (s, /^) ; ^oi^' '-^oi^ functions of order 
Oi{e, jj), depending on u, hio, /oi, g'oi and on the normal form transformation; Q^^^ : — > 
R^ is an ^-dimensional vector function, related to uj{Y) by 

N in 

) (Y; e,^l)^oo{Y) + J2Yl ^3,m-AYy l^"^'' , 

m=0 j=Q 
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being flj^rn-j{y) known vector functions and being Qqo = 0. The meaning of (5) is that 
the nomalized equations take the following form: the X-variation is provided by a modified 
frequency fi^^-* plus higher order terms in the normalization order or in the non-resonance 
order K; the variation of the normal form variable Y is constant, beside higher order terms in 
N and K. 

The coordinate transformation E^^^ results from the composition of a transformation Sc^^ 
acting on the conservative part and a change of coordinates S^^^ acting on the dissipative part, 
say: 

(r,X,t) = 5f)oSf)(|/,x,i). (6) 

Setting the intermediate variables as {y,x,t) = Ei^\y,x,t), the transformation E^^^ is im- 
plicitly defined through a sequence of generating functions i/jjo — i/^joiy, x, t), j = 1, N, 
as 

N 

X = x + ^ipjo^y{y,x,t)e^ 

N 

while E^^^ is defined by introducing suitable functions ajm, Pjm, j, £ '^+, as 

N m 

m=0 j=0 

N m 

m=0 j=0 

with the properties that aio{y,x,t) = /3io{y,x,t) = for any ?' > and {aij{y, x,t))s:,t = 
{/3ij{y,x,t))x,t = 0. The proof of the following lemma provides an algorithm to compute 
explicitly the vector functions ijjjo, ajm, /3jm, together with a suitable drift function rj, which 
allows to achieve the desired normal form. 

Normal Form Lemma. Consider the real analytic vector field (2) defined on A x T^+^ with 
complex extension in C,.„(A) x Cs(,(T^"'"^) for some Tq, Sq > 0. Let D C A be such that for 
any y e D, the frequency u = u{y) satisfies (3) for some K e Z+, a > 0. Then, there exist 
Eo, Ho > depending on ro, sq, K and the norms of u, h, f, g and there exists rj — rj{y, x, t), 
such that for (e, < (eq, fXo), one can find a coordinate transformation close to the identity, 
say S(^) : A x T^+^ ^ x T^+^ for a suitable normalization order N e Z+, which brings 



= x + ip!f\y,x,t) 

= y + i^i''\y,x,t), (7) 
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(2) into a normal form of order N as in (5). For fo < Tq, the normalized frequency Q^^^ in (5) 
is bounded by 

where A = max(£, \ and C^i is a suitable positive constant depending on tq, N and on the 
norms of u, h, f, g, respectively. Let Rq < tq, Sq < sq; with reference to the normal form 
equations (5), the following estimate holds: 

\\GN+l\\lh,So + \\Go-^\\ro,So + ^\\hii^J\Ro,So + l/^llboi^lko.So < CyX^'^^ , (9) 

for some positive constant Cy depending on ro, sq, N, K and the norms of u, h, f, g. 
Choosing^ N — [Kro/\ logAj] for some tq > 0, one gets that (9) becomes 

\\Gn+i\\r^,So + \\Go^\\ro,So + ^\\hi(^x\\Ro,So + l/^llboi^lUo.So < CyAe"'^'"° . 
Finally, denoting by Uy the projection on the y-coordinate, one has 

||n,(5f)o5f))-/4<C,A, (10) 
for some constant Cp > depending on tq, Sq, N and on the norm of u, h, f and g. 

Remark. As an outcome of the proof of the Normal Form Lemma, the drift i] will depend 
only on the normal form variable Y, i.e. r] = r]{Y). The explicit form of i] depends upon the 
functions hio, /oi, goi appearing in (2); we remark that the value of rj is determined by the 
requirement that Y is constant up to the normalization order A^, say Y — Yo + On+i{s, jj)', on 
the other hand, Yq depends on dq, which is chosen so that the frequency u — u{y) satisfies (3). 
The fact that r] is linked to the form of the vector field and to the frequency of motion appears 
also in KAM proofs for dissipative (or conformally symplectic) systems (compare with [13], 
[11]). 



Before giving the proof of the Lemma, we provide the statement of the main result, namely 
the confinement of the y variables for exponential times, which will be obtained through the 
Normal Form Lemma under the non-resonance condition (3). 

Theorem. Consider the vector field (2) defined on Ax T^+^ , and letD C A be such that for any 

y E D the frequency cu = uj{y) satisfies (3). Assume there exists sq, hq such that for {e, < 
(£o, A*o) the Normal Form Lemma holds. Then, there exist positive parameters po,To > 0, such 

^The choice of N can be performed as follows. The relation = e"^"^" impUes AT log A = —Ktq, namely 
iV = [Kto/\ log A|], where [•] denotes the integer part. 
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that for every solution {y{t),x{t)) at time t > with initial position (y(0), x(0)) e D x 
one has for X = max(£, |//|); 

\\y{t)-yiO)\\<poX for t < Cte''-\ 

for some positive constant Ct, where po and Ct depend on ro, sq, N, K, eq, /xq and on the 
norms of u, h, f, g. 

Remark. Notice that we obtain the standard formulation of the stability time in terms of an 
exponential estimate in the inverse of the small parameters, by adopting a proper choice of 
Ktq, say Kto < (^)'^ for some constant c > 0; in this case one has that the stability time 

estimate is t < Ct ^ for Aq = max(£o, Mo)- 



4 Proof of the Normal Form Lemma and of the Theorem 

In this section we first outline the general scheme of the proof and then we provide the com- 
plete proof of the Normal Form Lemma, followed by that of the main Theorem. For easiness 
of readability technical Lemmas appear in the Appendixes. We start by implementing a coor- 
dinate change of variables of the form (6). In particular, we define the intermediate variables 
(y, X, t) e R^x T^+^, which provide the transformation (7) in order that the following conser- 
vative normal form is obtained: 

N N 

i - ^?'\y-,e) + pj:Fjl'''\y,x,t)e^ + Y.Fji''''>^\y^ 

3=1 i=i 

+ ^h^o!y{y{y, S:, t),x{y, x, t),t) + pfoi^{y{y, x, t),x{y, x, t),t) + Ojv+i(£, p) 
y = pF^^'^\y,x,t) - fi{riN-i,i{y^S;,t)e^~^ + ... + rioN{y,x,t)p^~^) 

- ^Ko^xiyiy^ ^> t),x{y, X, t),t)+ pQoi^ {y{y, x, t),x{y, x, t),t) 

N 

+ E Fjo'''^''\y^ ^' t)e^ + ON+i{e, p) , (11) 

for suitable functions Fj|'^\ F^^'^K Fj}'''>''\ and being = u;{y) + 

J2f=i ^jo{y)£^, where ftjo{y) can be explicitly determined. We denote the inversion of (7) as 

X = x{y,x,t) = X + T^'''^\y,x,t) 

y = y{y,x,t) = y + T^y'''\y,x,t), (12) 

which provides the transformation Ei^\ We will see that the generating function ipjoiy, x, t) 
in (7) at the generic order j is the solution of an equation of the following form, defined in 
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terms of the intermediate set of variables: 

<^{y) '0io,x(^, X, t) + ^jo,t(^, X, t) + U^'-^\y, x,t) = , 

for a suitable known function L^^'-^^ (y, x, t) with zero average over {x, t). The above equation 
can be solved provided to = u{y) satisfies a non-resonance condition of the form 

k ■ u{y) + m ^0 for all e , m e Z , \k\ + \m\ < K , 

which is guaranteed by (3) on a suitable domain. 

After the implementation of the conservative transformation, to achieve the normal form (7) 
we construct a change of coordinates {Y, X, t) = S^^-* {y, x, t) defined as in (8), which allows 
to obtain the normal form (5). We will see that the functions /3jm must satisfy an equation of 
the form 

uj{Y)l3,mAy^ X, t) + X, t) + Nf^{Y, X, t) - X, t) = O , (13) 

for some known function NyrniX-, X, t); therefore, equation (13) can be solved provided that 
the drift components r}jrn{Y) are chosen as the averages Nj^ of N^^'- 

The normal form equation (13) can be solved provided that the frequency satisfies the non- 
resonance condition k ■ U!{Y) + m 7^ for all /c G Z^, m e Z, \k\ + \m\ < K, which is 
guaranteed by (3) on a suitable domain. Once /3jm is determined, we can proceed to compute 
ajjn by solving a normal form equation again of the form (13), but having zero average. 

Proof of the Normal Form Lemma. We prove by induction on the normal form order that 
we can determine the transformations (7) and (8) so to obtain the normal form (5). We start by 
constructing the first order normal form through the implementation of the conservative and 
then of the dissipative transformation; in a similar way we construct the transformations for 
the order N. Being the proof quite long, for sake of exposition we split it into separate steps. 

Step 1: Conservative transformation for iV = 1. 

Let us start with the conservative transformation for N = 1, namely we implement the first 
order change of variables 

X = X + e^ioAy,x,t) 

y = y + eil^ioAy.x.t) , (14) 

where ^10 = V'io(^, x, t) must be determined. Let fo < ro, 5q < sq, sq = sq — So; then (14) 
can be inverted as 

x = x + er^'''^\y,x,t) 
y = y + er^^''\y,x,t), 
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for suitable functions F*^^'^) and F^^'^^, provided the following condition is satisfied (see Ap 
pendix A): 

""'-'10,1/ ||r(),so 6 

Using (14) and (2), we compute the time derivatives of x, y as 



rO£l|^io,J|r„,.oe'^°V'<l- (15) 



X = uj{y) + eujy{y)tlJio,x{y,x,t) + ehio,y{y,x,t) + ehio,y{y) + lifoi{y,x,t) 
+ euj{y)ipio,yx{y: X, t) + eipio -ytiy, x, t) + O^ie, n) 

y = -£/lio,a;(^,5,t) +/i(fi'0l(y,5,^) -^Ol(y,^,^)) -£w(y)^/'lO,a:a:(j',^,^) 

- £ipw,xt{V-: ^) + 02{e-, A*) • (16) 

The conservative normal form is obtained by imposing that V'io(y, x, t) satisfies the following 
normal form equations: 

i^y{y)i^io,x{y^ + i^io^ytiy, X, t) + uj{y)'ilJio,yxiy, X, t) + hf^yiy, x,t) = 

'^(y)^10,a;a;(^,^,i) + V'l0,xt(y,5,^) + = 0. (17) 

As a consequence, setting 

n^;\y;s)=iu{y) + ehio,y{y), (18) 

equations (16) become 

X = n^^\y;e)+eh^o^y{y,x,t)+ iifoi{y,x,t) + 02ie,iJ,) 

y = -eh^^^{y,x,t) + ii(^goi{y,x,t) -r]oi{y,x,t)^ +02{e,fi) , (19) 

which are recognized as being of the form (11). We remark that equations (17) are obtained 
taking, respectively, the derivatives with respect to y and x of 

w(y)V^io,x(y, X, t) + ^io,t(y, X, t) + hf^iy, x,t) = . 
Expanding hff^{y, x, t) in Fourier series as 

hf,''{y, x,t)= j2 ko,km{yy^'"'^"''^ , 

(fc,m)eZ^+Mfe|+|m|<i<: 

where hw,km{y) denote the Fourier coefficients, the solution for V'io(y, x, t) is given by: 

(jj(y) ■ k + m 
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To control the small divisors appearing in the previous expression, let us invert the second in 
(14) as y = y + eR^^\y, x, t) for a suitable function R^^^ — R^^\y, x, t), provided that for 
Tq < To one has (compare with Appendix A) 

70£||^io,x||r-o,.o:^^ < 1 ■ (20) 

To - To 

Then, we have that 

\uj{y) -k + ml > a - eK\\R'^'^^y^^sJ\ujy\\ro > ^ , 
if (compare with Appendix A) 

a 

' - 2K\m\,,M\r, ■ ^^^^ 

Step 2: Dissipadve transformation tor N — 1. 

We proceed to reduce to normal form the dissipative part through a first-order transformation 
of coordinates A^^\ which we write in components as 

X = X + aoi{y,x,t)iJ, 

Y = y + /3oi{y,x,t)n, (22) 

for some functions aoi and /3oi to be determined as follows. We premise that equations (22) 
can be inverted as 

y = Y + A^y^'\Y,X,t)ii 

for suitable functions A^^'^^ and A^^'^^ provided the following conditions are satisfied (see 
Appendix A): 

70|;u| ||aoi||fo,5o e^'°So^ < 1 
70|//| (||/3oi||fo,5o + l/^l 11/^01 ||aioi||fo,so) 5- < 1 , (23) 

where = So/2, Rq < tq and being || A^'^'^) \\ro,So < ||aoi||fo,5o with Sq < Sq - Sq. Up to the 
second order, the inversion of (22) provides 

X = X - aoi{Y,X,t)fx + 02{fi) 

y = Y-MY,X,t)f, + 02{fJ,). (24) 
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Taking the derivative of (24) and using (19), we express X,Y asa function of X, Y as 

X = uj{Y)-Uy{Y)f3oiiY,X,t)fi + eho,y{Y) 

+ eh>,'^y{Y, X, t) + /x/oi(F, X, t) + u{Y)aoi,.{Y, X, 

+ aoi,t(^,X,t)/i + 02(e,//) 

y = -eh>,^,{Y,X,t)+fx(^go,{Y,X,t)-Voi{Y,X,t)) 

+ a;(r)/3oi,.(r, X, t)fi + /3oi,t(r, X, t)fi + 02{e, fi) . 

The dissipative normal form is obtained imposing that ccqi, Poi and 7701 satisfy the following 
equations: 

u;iY)aoi,,iY,X,t) + aoiAY,X,t)-Uy{Y)f3oi{Y,X,t) + f^,''{YX,^ = 
u{Y)PoiAY,X,t)+My,X,t) + ~gi''iY,X,t)+goi(Y)-rjoi{Y) = 0, (25) 

where we have split /oi into the sum of its average and of the oscillatory part, namely 

/oi(y, X, t) = /oi(y) + fi^'iY, X, t) + f>''{Y, X, t) 

and similarly for ^01: 9oi(Y, X, t) = goi{Y) + g^^{Y, X, t) + ^01^ (F, X, i). From the first of 
(25) we see that the average of /3oi is zero and we can assume that also the average of aoi 
is zero. Then, equations (25) can be solved, provided that in the second equation the term 
7701 (y, X, t) is chosen so that 

r^oi{Y,X,t)=rioi{Y)=MiY). 

The final normal form can be written as 

X = u{Y) + eh^,y{Y) + /x/oi(y) + eh>^^{Y, X, t) + /x/o>^(y, X, t) + F,{Y, X, t) 
Y = -sh>,%{Y,X,t)+i,g>f{Y,X,t) + G2{Y,X,t), 

where F2, G2 are O2 {e, n). These equations are recognized to be of the form (5) with il^^^ (Y) — 
cu{Y)+sho{Y)+MY). 

The solutions of (25) involves small divisors, which can be bounded as follows: 

\uj{Y) • A; + m| > I - |/i|ir||/3oi||r-o,Soll^3/llr-o > ^ , 
provided that (see Appendix A) 

1^1 < AT^WR 11" II II ■ (^6) 
llPoi||fo,soll'^j/llro 
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Step 3: Conservative transformation for the order N. 

Assuming that the Lemma holds up to the order N — 1, we prove it for the order N, starting 
from the change of variables (7), that we invert as 



N 



X = x + J2^f^iy^^^ty -^NoAy,^,t>'' + ON+i{e) = x + r^^'''\y,x,t) 

AT 

y = y + J2^'g{y,x,t)e^ + ^IJmAy,S:,t)e'' + OM+iie)=y + r^y''^^ (27) 

(x) (v) 

where Tjq, F^q are known, since they depend on the known functions ^lo, ipN-ifl, while 
P(x,jv)^ p(y,iv) j^^^g ijggj^ introduced as in (12). Choosing fo < tq, Sq < Sq, the inversion is 
possible provided that (see Appendix A) 

70||^f)||.-o,.oe''°^o' < 1, (28) 
being ip^^'> — Ylj^=i i^jo^^- For short, let us write the equations (27) as 

x = x{y,x,t), y = y{y,x,t). 
Inserting (27) in (2) and expanding in Taylor series, one has 

X = Lj{y)+Uy{y)ijNoAy,S:,t)e'' + F^'''''\y,x,t) 

+ ^h^o!yiy{y^ ^> t),x{y, X, t),t)+ nfm^{y{y, x, t),x{y, x, t),t) + On+i{^, jj) 

y = G'(°'^)(y,£,t)-^(^^_i^i(y,5,t)£^-i + ...+r7o,jv(y,5,t)/i^"') 

- ^h^o!xiyiy^ ^> t),x{y, X, t),t)+ i^goi^{y{y, x, t),x{y, x, t),t) + ON+i{e, n) , 

where F^^'^), are known functions; contains terms of order e, e^, , ji, jie, 

fj,e^~^, while contains all powers ii"^ with 1 < j+m < N. Next step is to compute 

X, y as a function of x, y. Taking into account (7) and that by the inductive hypothesis ip^^\ 
^(•^-1) make the equations in normal form up to the order A'^ — 1, we obtain 



+ Yl ^jl^^\yy + ^y{y)i^m,x{y^ + ^{y)i^No,yx{y, x, t) + ipm,yt{y, x, t) 



N 

X — Uj[~] ' 

i=i 

N-l N 

t3 



i=i i=i 
+ ^Ko^yiyiy^ ^> t),x{y, X, t),t) + nf^^^{y{y, x, t),x{y, x, t),t) + Ojv+i(£, A^) , 
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where F(i'^)(y,x,i) is a known function that has been decomposed as 

i=i 3=1 

N-l 

In a similar way one obtains 

y = (y, X, t) + i^Pj^f'-^^ {y, x, t) - uj{y)ipm,xx{y, x, t) - ipm,xt{y, x, t)j 

N 

- /i(77^r-i,i(y,x,i)£^-^ + ... + 77o,iv(y,x,i)/.^-^) + J]i^o''^^^^^ 

- ^Ko^xiviy^ ^> t),x{y, X, t),t)+ idg^i^{y{y, x, t),x{y, x, t),t) + Ojv+i(£, 1^) , 



for known functions ^(3,a^)^ ^j^e latter having zero average. The conservative normal 

form is obtained by imposing that V'ivo solves the following normal form equations 

u}y{y)ipNo,x{y^^^t) +^{y)'^No,yx{y,S;,t) +ilJm,yt{y,S;,t) + F^^f'-^\y,x,t) ^ 

i^{y)'ipm,xx{y^ t) + i^N(ixt{y^ t) - FNf'~^\y^ ^> ^) = o (29) 

where as before we have split the known function -F^o^'' into -F^o^'' = -^jvo^'' + -^jvo^'' well 
as -F^o^'' into F^^^'^^ = F^q^'~'^^ + F^o^'^^-*. Note that in this setting the A^-th order contri- 
bution to the shifted frequency vector is given by QNo{y) = F^^q^^ (y)- Due to the Hamiltonian 
character which occurs for /x = 0, there exists a function M^^'-^^ — M^^'-^\y,x,t) with 
zero average, such that 

Qy —^m {y,x,i) , — -r^o {y,x,i) , 

so that (29) are equivalent to solve the equation 

uj{y)i^m,x{y^ ^> t) + i^m,t{y^ + = . (30) 

The solution of (30) provides the function ipNoiy, which produces the conservative nor- 
mal form: 

N-l N 

3=1 3=1 
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+ ^h^o!y{y{y^ 3:, t),x{y, x, t),t)+ f^foi^{y{y, x, t),x{y, x, t),t) + ON+i{e, n) 
~ ^ ^F^'''^\y,x,t) - i^{r]M-iAy,S:,t)^'''^ + ■■■ + Vo,Niy^S:,t)i^''-') 
- ^h^o!x(y(y^ ^> t),x(y, X, t),t)+ HQoi^ (y(y, x, t),x(y, x, t),t) 



N 



+ 5Z^o''^'^''ny,^,^y+ 0^+1 (31) 

where 

e) ^ u{y) + F\l^''\y)e^ . (32) 

Notice that from the definition (32) one obtains an estimate Uke ||Jlc^^ ~ ^|| < C'cS for a 
suitable constant Cc- The solution of (30) involves small divisors of the form uj{y) • k + m, for 
k G Z^, m G Z with + |m| < /T, as it was for the case N — 1; non-resonance is guaranteed 
provided that (see Appendix A) 

£ < ...M^.^n," n „ , (33) 



where R^^'^ is the function such that y — y + eR''^\y, x, t) and fg < fo. 
Step 4: Dissipadve transformation for the order N. 

As for the dissipative part, we consider the transformation (8), that we invert as 

N N~k N-1 

fc=0 3=0 k=0 

N N-k N-1 

y = ^ + J2Y1 ^' - Yl ^k,N-k{Y, X, t)£V^"' + ON+i{e, ii) , (34) 

fc=o i=o fc=o 

for suitable known functions akj{Y, X, t), bkj(Y, X, t) with a^o = bko — for A; = 0, N, 
provided the following conditions are satisfied (see Appendix A): 

70(||^(^)|ko,.-o + ll/9f^l|.o,.-oll«^^^l|.o,.-o) < 1, (35) 

To — rio 

where = so/2, -Ro < ^o- For short let us denote (34) as 

x = x{Y,X,t), y = y{Y,X,t), 
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while we express the original variables in terms of {Y, X, t) through 

X = x(y(Y,X,t),x{Y,X,t),t)=X + ^^'''^\Y,X,t) 
y = y{y{Y,X,t),x(Y,X,t),t)=Y + ^^y'''\Y,X,t), 

for suitable functions $(^'^), $(2^'^), which are /x). We need to determine the unknown 
functions Q!o,iv, aN-1,1, Po,n, ■, Pn-i,i, VN-1,1, r]o,N as follows. Starting from (31), we 

compute X, y in terms of X, Y and we express X, Y in terms of X, Y, using (8), (31), (34); 
by the inductive hypothesis a^j, f3kj, i]j with 0<A; + j<A^ — 1 are determined so that the 
equations of motion are in normal form up to the order e'^n^ with 0<k+j<N — 1. This 
leads to the following equations: 

(N-l \ N-1 

J2 f3j,N-j{Y, X, t)e^^x''-^ + u{Y) J2 o^j,N-Uy, X, i)^^^-^' 

N-l 

+ «i,iv-,,t(>^, X, t)£>^-^- + X, t) 

3=0 

N 

+ E F\l'''^^''\m X, t).x{Y, X, t),t)e^ 

+ eh>,l{y{Y, X, t),x{Y, X, t),t) + /x/o>i^(2/(F, X, i), x(y, X, t), t) + 0^+i(£, /x) 

AT-l 

N~l 

3=0 

N 

i=i 

- eh>,'^MY, X, t),x{Y, X, t),t) + iig>^^{y{Y, X, t),x{Y, X, t),t) + O^+iie, /i) , 

where F^^'^^Y, X,t), F^^'^\Y,X,t) axe 0(/i^-i, e/x^-^, e^'V, ^^"')- Let us decom- 
pose as X, t) = F(^'^)(F) + X, t). The normal form at the order 
N is achieved imposing that akj, hj, Vj satisfy the following normal form equations 

-ujy{Y)Pj,N-j{Y, X, t) + u{Y)aj,N-jAY, X, t) + a,,N-j,tiY^ X, t) 

+F^i:f\Y.x^t) = ^ 

u{Y)P,,^_UY, X, t) + P,,N-,,t{y^ X, t) + F\%''_:f\Y, X, t) 

+fI%''-]{Y) - V3,N-j{Y, X, i) = (36) 



18 



for < j < N — 1, where we have used the expansion 

N-l N-1 
fc=0 k=0 

and we split Fjf^^l = Fun^'^^^ + ^kN-'k"^^- non-resonance condition for a'(y) reads 
as (see Appendix A): 

K\\l3^''^\\ro,soH\\ro<l, (37) 

where, for short, we have written Y = y + x, t; £, //). From the second of (36), it is 

so that the second of (36) can be solved to determine Pj,N-j, while from the first of (36) we 
obtain aj^^^j. Identifying the final normalized frequency with n^f^\Y ; e , fj,) = oj{Y) + 
Em=oEiLo^-,rn-i(^)£-'Ai"'"^' , whcrc ^00 = 0, 1^10 = uj, Qjo = F^^'^\y) and ^j,m-j = 
-P^^ti-j(^)' resulting normal form is given by 

N-l 

X = Of)(F;£,/.) + Fjv+i(y,X,i) + 5^i^>^f^)(F,X,i)£>^-^- 

j=0 

+ sh>,l{y{Y, X, t),x{Y, X, t),t)+ /x/o>i^(2/(F, X, t),x{Y, X, t),t) 

N 



+ E Fll''''^''\yO'. X, t),x{Y, X, t\t)e^ 



N-l 



Y = G^+^{Y,X,t) + Y,F'i%''--''\y,X,t)s^l^^~' 

- ehX^MY^ X, t),x{Y, X, t),t) + f^g>f{y{Y, X, t),x{Y, X, t),t) 

N 

+ E F^o'''''''\y{Y, X, t),xiY, X, t),t)e^ , (38) 
where Fjv+i, Gat+i are Ojv+i(£, //). For short, we define 

N-l N 

F>^{Y, X,t) ^ Yl ^' t)e^^'^~' + E X, t),x{Y, X, t),t)e^ 

j=0 j=l 
N-l N 

G>,''{Y, X,t) ^ ^' + E Fio'''''''\m X, t),x{Y, X, t),t)e^ , 

j=0 j=l 
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which makes (38) of the form (5). The estimate (10) follows from the fact that (7) is close to 
the identity up to terms of order e, while (8) is close to the identity up to terms of order /x. 
Notice that the smallness requirements on e, fi, say s < Eq, < /iq, are needed to ensure that 
the non-resonance condition (see (21), (26), (33), (37)) is satisfied and that the transformations 
(7), (8) can be inverted (see (15), (20), (23), (28), (35)). 

The original variables can be expressed in terms of the intermediate variables by means of 
(27), provided (28) holds with 

II J- l|ro,so — II rj/ \\ro,so ) 

being so = so — Sq. Moreover, we have 

||J- ||ro,so — lira; \\ro,so ^ WYxx Iko.soll^ lko,so ■ 

For the same reason, equations (34) are invertible for £, |//| sufficiently small, since the Jaco- 
bian of the transformation is close to the identity; we can write the inverse as 

X = X + A^'''^\Y,X,t) 
y = Y + A^y^^\Y,X,t) 

for suitable, bounded functions A^^'^^, A^^'-^) of order Oi{e, ji). We finally obtain 

X = x + A(^'^)(r,x,t) + r(^'^)(y + A(^'^)(r,x,t),x + A(^'^)(r,x,t),t) 

= x + $(^'^)(y,x,t) 

y = y + a(^'^) (y, X, t) + r(2''^) (y + a^^-^) (y, x,t),x^ A^^-^) (y, x, t) , t) 

= Y + ^^y'^\Y,x,t) . 

Recalling Lemma B.l of Appendix B, for Tq > and 1 < j < ^, we have 

\\hi^xyj\\Ro,So ^ C'a||^10,a;j/jRo,S'o+roe~^^'*"^^'^ 

\\h>^ II <^ r \\h II ^-{K+1)to 

\\'^10,xxj\\Ro,So S '~^a\\n'10,xXj\\Ro,So+To<^ , 

where Ca is a positive real constant. Setting 

C(<s>,y) = sup ||<l>5.^'^^||i?„,5o , C^<s,^x)= sup \\<^f'^^\\R^,So , 

i<j<e i<j<i 

C(h,y) = sup ||/ilO,xyJ|flo,So ' C{h^x)= SUp ||/ilO,a;a;Ji?o,So ' 

i<j<i i<j<e 
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one obtains 



\h>,%{Y + <^(y'''\Y,X,t),X + <!>(^''^\Y,X,t)MRo,So < \\h^^^^^ 



+ e( sup ||/iio,c.%||Ro,5o sup \\^f'^^\\R^,So+ sup \\hw,xxj\\Ro,So SUp '^^ ||ro,5o ) 

\T<-i<f l<j<£ T<i<l> T<i<f / 



i<j<e 



l<j<E 



i<j<e 



having defined Ch = CaC ||/iio,x||Ro,5o+ro + + • Similarly we 

obtain 

\\g>f{Y + $(^'^)(y, x,t),x + $(-'^)(y, X, i), i)|Uo,5o < C.e-^^'' , 

where 



sup ||S'oi,2/3l|i?o,So ' 



C{g,x) = sup ||5'oi,a;J|flo,So ' 



Analogously we find 

||G>i^(F + $(^'^)(F,X,t),X + $(^'^)(F,X,t),t)||«„,So < ACce-^^" 

where 



C(G,2/) = sup ll/^o ^GqI^. ||ro,So 



C, 



G 



C(G,a;) = sup ||/^o ^^OM, II Ro,So 



Let us bound Gn+i in (38) as 

for a suitable constant Cg- From the second of (38) we obtain: 



|<^Af+l||i?o,So + ^II^1o!^IIro,S'o + l/^IIIS'oi' 



I -Roi'S'o 



Kto 



||<-^01 IIRo.So ^ <-^gA 



Choosing N as 



log A I 



we have (see (9)) 

||GAr+l||ilo,5o +s||^l{M:l|iio,So + l/^lllfoi^llflo,So + ll^^Ol^llflo.So ^ CyX^^^ , 



(39) 
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with Cy = Cg -\- Ch-\- Cg -\- Co- This concludes the proof of the Lemma. □ 

Proof of the theorem. The distance between the solution at time t > 0, say y{t), and the initial 
condition |/(0) can be bounded as 

\\y{t) - ym < im - rm + \\Y{t) - rm + \\y{o) - ym . 

By the estimate (10) of the Normal Form Lemma one has 

\\y{t) - Y{t) \\<C,X, \\y{0) - y (0) II < C^A . 
By the second of (5) and by (9), one has 

\\Y(t) - y(0)|| < r {\\G^^r\U,So + e\\h>,%\U,So + l/^llbo>i^ll«o,5o + \\G^i''\\ro,So) ds 
Jo 

Let r2 > be such that 

t Cy\^ < r2 , 

which is satisfied as far as 

t < -If A-^ < e^-° , 

Oy Oy 

where (39) has been used. Finally, setting po = 2Cp + r2, we obtain 

||yW-y(0)||<poA for i<ae^^°, 
having defined Ct = r2/CY. O 

5 An application of the normal form and of the stability es- 
timates 

In order to test the accuracy of our results, we implement the Normal Form Lemma and we 
derive the stability estimates on a specific example. To this end, let £ = 1 and let us consider 
the differential system: 

X — y — /J. (sin(a; — t) + sin(a;)) , 

y — —e{sm{x — t) + sm{x))—ijL{y — r]). (40) 

We remark that for = the system (40) is associated to the Hamiltonian function in the 
extended phase space 

H{y, T,x,t) = — — s {cos{x — t) + cos{x)) + T , 
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where the unperturbed frequency of the motion is given by uj{y) = y and being T conjugated 
to time. We provide details for the computation of the second order normal form associated to 
(40) (see Section 5.1). A comparison with a numerical integration is performed in Section 5.2. 
Stability estimates according to the Theorem of Section 3 are computed in Section 5.3. A 
slightly different example with oscillating energy is analyzed in Section 5.4. 



5.1 Normal form 

The second order normal form can be computed as follows ^. At first order we identify the 
non-zero average contributions by hio,y{y) = (see (18)). The conservative normal form 
equations become (see (17)): 

y V'lO,a;a;(y,^,^) + ^10,a;t(^,5,^) +Sin(x-t) +Sin(£) = 0, 

from which we get 

, , _ , sm(x — t) sm(x) 

y - 1 y 

The second order conservative normal form is obtained by computing the generating function 
11^20 as the solution of the equations: 

cos(2x — 2t) 

V'20,x(y, X, t) + ^20,ytiy, X,t)+y ^20,yxiy, X, t) + _ 

(l-2y)cos(2x-t) cos(t)(l-2y) cos(25) _ 
+ 2(^-1)2^2 + 2(y-l)V W~ ~ ' 

sin(25; — t) sin(25; — 2t) sin(25;) 

y^2o,xx{y,^,t) + ^2o,xt{y,^,t)+ ^y_i)y + 2{y-iy " 

which provides the function V^2o as 

sin(2a; — t) sin(25; — 2t) sin(t) sin(2£) 



i^2o{y,x,t) 



2{y - l)y{2y - 1) 8{y - If 2{y - l)y 8^ 



while the second order term of the frequency shift is given by i^2o{y) = (— + Sy^ — 3y + 
l)/(2(y - l)^y^), so that ^l^c\y) = oj(y) + Jl2o(y)£^ (see (32)). At this stage, we succeeded 
in normalizing the symplectic contributions and in getting the conservative normal form to the 



^Notice that here we first implement the conservative transformation to the second order and then we deter- 
mine the dissipative transformation, which provides the second order normal form. 
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second order in the intermediate variables as 

X = -f::-^-^,:^^e^ - i^{sm{x - t) + sm{x)) + Osie) 

y = -Ky -v) + Osie) . 

The first order of the dissipative normal form, expressed in terms of the new variables, provides 
the equations 

yaoi,.(F, X, t) + aoi,t{Y, X, t) - /3oi{Y, X, t) - sm{X -t)- sm{X) = 
Yl3oUY,X,t)+l3oi,t(Y,X,t)+Y-rioi(Y,X,t) = 

(see (25)). From the second equation we get 

r}o,{Y,X,t) = Y , (3oi{Y,X,t) = 0, 

while from the first equation we obtain 

^ fV Y +\ C0S(-^ - t) COS(X) 

In a similar way, the second order dissipative normal form provides the equations 

Fan,.(r, X, t) + an,i(F, X, t) - f3uiY, X, t) + ^^^y'^^^'y^^ = 0, 

YPnAY, X, t) + PuAY, X, t) + Vu{Y, X, t) + ^^P^ - ^'"^ '2(F -1)^ " ' 
cos(2X - 2t) ^ (1 - 2F) cos(t) cos(2X) _ ^ 



2(Y-1) 2{Y-1)Y 2Y 

-1 

2(y-i)y ' 2{Y-i) 



Yo^AY, X, t) + a^AY, X. t) - MY, X, t) + '^^ -JJ'lf J " + ^^^^ + 



(1 - 2Y) cos(i) cos(2X) 
2{Y-l)Y ^ 2Y " ' 
YP02AY, X, t) + ^^2,t{Y, X, t) + rio2(Y, X,t)^0. (41) 

From the second and fourth of (41) we conclude that 

V02{Y) = 

2Y -1 

Vu{Y) - 



f3ii{Y,X,t) = 



2{Y - 1)Y 

sm{2X-t) sm{2X-2t) (1 - 2y) sin(i) sin(2X) 



2(1 -y)y 4(y-i)2 2(y-i)y 4y^ 
MY,x,t) = 0, 
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while from the first and third of (41) we obtain 



_ cos(2X-t) , cos(2X-2t) , 
«iiir,A,tj - 2y(2y2_3y + i) + 8(y-l)3 + 

(-2y^ + + 3y - 2) cos(t) cos(2X) 

2{Y - i)2y2 ^ 

sin(2X - t) sin(2X - 2t) {2Y - 1) sin(t) sin(2X) 

ao2{Y,x,t} - 2(y - i)Y ^ 4(y - 1)^ + 2(y - i)y ^ w " 

In conclusion, the second order normal form associated to (40) is given by 

• {1-3Y + 3y2 - 2y3) 2 (1 - 2y) 2 ^ , , 

y = o + 03(£,/^), 

where 77 = vO^) takes the following expression: 

This concludes the computation for the second order normal form. In a similar way one can 
continue to higher orders and in fact we computed up to the fifth order were we stopped, since 
i) we already reached exponential estimates for non trivial parameter values (see Section 5.3) 
and ii) the structure of the terms to be determined becomes too complex to allow for higher 
order computations using just a general purpose algebraic manipulator (Mathematica 7); nev- 
ertheless, we believe that higher orders can be obtained by implementing a specific algebraic 
manipulator in C or Fortran languages. In order to understand the degree of complexity of the 
computation, let us denote by Sc^^ o S^^^ respectively, the conservative normal 

form at order A^, the dissipative transformation and the overall normal form. At any order N 
the algebraic manipulator has to deal with Poisson series ([22]) of the form 



' Qikiy) 



E 



where [/ is a sublattice of ajk are complex coefficients, bjk, Cjk G Z+ with bjk + Cj^ = N 
and Pji^;, Qjk are polynomials in the actions. Let us denote by deg{Pjk) the degree of the 
polynomial Pjk, which contains only positive (or zero) powers in the action (and similarly for 
Qjk)- For each order N between 1 and 5 the numbers of Fourier terms as well as the degree of 
the polynomials Pjk and Qjk are provided in Table 1 . We remark that the main limitation in the 
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Table 1 : Number of Fourier terms and the degree of the polynomials of the conservative and 
dissipative transformations as a function of the order. 







N^2 


A^ = 3 


A^ = 4 




# Fourier terms for Sc^^ 


9 


26 


130 


524 


1888 


# Fourier terms for S^^^ 


6 


25 


201 


846 


6829 


# Fourier terms for Sc^-* o S^^-* 


11 


53 


461 


2875 


5004 


deg{Pjk) for 





4 


6 


8 


10 




1 


3 


5 


7 


9 


deg{P,,) for Ef^oEf^ 


2 


5 


6 


18 


29 


degiQjk) forSi^^ 


2 


1 


1 


1 


2 


degiQjk) fovEf^ 


1 


1 


1 


1 


1 


degiQjk) for si^^ o Ef^ 


2 


2 


2 


11 


19 







present implementation of the normal form algorithm turned out to be the capability of dealing 
with the algebraic manipulation of fractions of polynomials of higher degree (29 at order 5)^. 



5.2 Comparison with a numerical integration 

To compare the results with a direct numerical integration we need to find the expression 
for r]{Y) in terms of the original variables (y, x, t). To this end, we determine ri{y) from the 
condition that still y = in terms of the original variables up to the normalization order N, 
thus obtaining 

Remark. By induction on the normalization order, one can easily prove that the normal form 
equation in new and old variables keeps the same functional form. For this reason also r] 
maintains the same form in old and new variables. 

^Note, that the number of tenns of a Poisson series expansion with rational coefficients in the actions strongly 
depends on the 'algebraic normal form' of them. Putting terms of the same Fourier mode under the same de- 
nominator will increase the exponent order in the denominators and reduce the number of Fourier terms. On the 
contrary, writing the sums apart will reduce the exponent order in the denominator, but increase the number of 
terms. The numbers given in Table 1 strongly depend on the choice of the form of the rational coefficients. 
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10-'° 
10-'^ 

10-14 



6 10' 1.2 lO* 1.8 ■lO'' 2.4-10'' 3 - 10* 

t 



Figure 1 : Comparison between the analytical and numerical solutions obtained using normal 
forms of 1st (upper curve), 3rd (central curve) and 5th (bottom curve) orders. The integration 
time is T = lO^n, the parameters are £ = 10^"^, /i = 10^^; the initial condition is set to 

ro = l(v/5+l),Xo = 0. 



For given initial conditions (Xq, Yq) we integrate the normal form equations (up to the nor- 
malization order) as 

X{t) = Qf\Yo;e,fi)t + Xo 
Yit) = Yo, 

where i^^^-* = Q.^J^\Y; e, ji) denotes the normalized frequency to the N-\h order. Then we 
back-transform to old variables and we define the relative error between the analytical and the 
numerical solution as 

^^^^^^ ^ 1 ((Xnum(t) - Xana(t)) + {VnnM - V-^nM) ^) ^''^ ^42) 
^ Valium 

where (xnumj Z/num) IS the state vector at time t obtained from a numerical integration of the 
original equations of motion and (a^ana, I/ana) represents the state vector at time t obtained 
from the normal form solution back-transformed in the original variables. The evolution in 
time of the error for the parameter values e = 10~^, /i = 10^^ and the initial conditions 
^0 = llv^S + 1), Xq = 0, is shown in Figure 1. We plot the value of err(t) versus time for 
an overall integration time of T = IO^'tt. The analytical solution was computed using the 1st, 
3rd and 5th order normal form. The numerical solution was obtained using a 4th order Runge- 
Kutta integration scheme with fixed step size St = 10^^. As expected, the difference between 
the numerical and analytical solutions decreases as the order of the normal form increases. 
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5.3 Exponential stability estimates 



In this section we present an application of the Theorem of Section 3 to the sample provided 
by the differential system (40). We first discuss the smallness conditions required for the pa- 
rameters (Section 4.3.1) and then we compute the stability estimates (Section 4.3.2). 

5.3.1 Bounds on the parameters 

The bounds on the parameters e and n are due to the smallness conditions imposed by the 
requirements to invert from original to intermediate variables, to invert from intermediate to 
new variables, to satisfy the non-resonance condition in the intermediate variables and the 
non-resonance condition in the new variables. With reference to the Appendix A, assuming 

= i (v^S + l), xq = 0, ro = 0.118, sq = 0.1, 5q = 0.05, K = 20, one finds fo = 0.113, 
sq = 0.05, fo = 0.056, i?o = 0.057, 5*0 = 0.025, a = 0.09 (the parameters are chosen so to 
optimize the result). Condition (28) requires that e < 1.2 ■ 10"'', while condition (33) imposes 
that e < 1.2 ■ lO^'' (no requirements are needed on n). Condition (35) is satisfied provided 
e < 1.2•10-^ 1/^1 < 2.0•10-^ while condition (37) requires that e < 3.0-10-3,r/ < 4.75-10-3. 
In conclusion we obtain that all conditions are satisfied provided that e < 1.2 • 10"^ and 

< 2.0- 10-^ 

5.3.2 Stability estimates 

The final step is to implement the estimates derived from the Theorem, keeping in mind that 
there are no Fourier modes of the form h^^ in the sample (40). Let us write the transformation 
of the second component from original to final variables as F = y + T^^\y, x, t). For given 
£ < ^0. Ia'-I < /^o and given Aq = max(£:o, A'-o). we calculate the constant Cp as ||7'*'^-'||ro,so/ -^o 
and we define ri = C^Aq. Taking r2 of the same order of magnitude as ri, we set Ct — r^jCy 
with Cy = IIGat+i ||i?o,So/ -^o^^^- For the variation in action space we set po — 2ri + r2Xo. We 
finally compute Tq from (39). In conclusion we obtain that 

\\y{t)-y{0)\\<po for any i < Tq = C^e^^" . 

The actual values of the parameters, as a function of the normalization order A^, are summa- 
rized in Table 2. The parameters Eq and /iq are taken from the estimates on the smallness of the 
parameters (see Section 5.3.1), while K was set equal to 20. 

Similar estimates can be obtained by fixing tq and calculating K accordingly. Since K also en- 
ters into the denominators of the non-resonance conditions (33) and (37), it will also influence 
the bounds on the smallness of the parameters. The results are shown in Table 3. Fixing tq — 2 
the choice of K depends on the order of normalization and on the bound on the smallness 
parameter Aq, leading to a slightly different Ktq compared to the values given in Table 2. The 
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N_ 2 3 4 5 

Vo 0.851 1.277 1.703 2.129 

\\Gn+i\\ro,So 1.966-10-^ 5.147-10-1° 1.320-10-12 2.053 - lO^^^ 

\\TN\\ro,so' 3.815-10-^ 3.819-10"^ 3.819-10"^ 3.820-10"^ 

Cp ' 1.908 1.909 1.909 1.910 

Cy 4.915 6.433-101 8.251 - 10^ 6.416 - 10^ 

Ct 3.881-10-1 2.968-10-2 2.314 - 10-^ 2.977-10-^ 

po 1.145 -lO--'^ 1.146-10-3 1.146-10-3 1.146 - 10-^ 

Tn 9.702 - 10^ 3.710 - 10^ 1.446 - lOi^ 9.302 - lOi^ 



Table 2: Stability results versus the normalization order N. For the definition of the constants 
see the text. The parameters are taken to optimize the stability time: Eq = 1.2 - 10-^, /xq = 
2.0 - 10-\ I/O = i(V5 + 1), ro = 0.118, Sq = 0.1, K = 20. Notice that in this table we fix K 
and we let Tq vary. 



N 


2 


3 




4 


5 


K 


8 


12 




17 


21 


\\Gn+i\\ro,So 


1.970 - 10-^ 


5.169-10- 


-10 


1.329 - 10-12 


2.072 - 10-1^ 


||ro,so 


3.815 - 10-^ 


3.819 - 10- 


■4 


3.819 - 10-^ 


3.820 - 10-^ 


Ao 


2.0 - 10-^ 


2.0 - 10-^ 




2.0 - 10-^ 


2.0 - 10-^ 


Cp 


1.908 


1.909 




1.909 


1.910 


Cy 


4.924 


6.461 - lOi 




8.307-102 


6.476 - 103 


Ct 


3.874 - 10-1 


2.955 - 10- 


-2 


2.299 - 10-3 


2.949 - 10-^ 


Po 


1.145 - 10-3 


1.146-10- 


-3 


1.146 - 10-3 


1.146 - 10-3 


To 


3.443 - 10^ 


7.828 - 10^ 




1.341 - 1012 


5.129 - lOi^ 



Table 3: Stability results versus the normalization order A^. For the definition of the constants 
see the text. The parameters are taken to optimize the stability time: Sq — 1.2 - 10-^, no — 
2.0 - 10-"^, yo = |(V^ + 1), ro = 0.118, sq = 0.1, tq = 2. Notice that in this table we fix tq 
and we let K vary. 
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0. 



1 2 3 4 5 6 

X 

Figure 2: Numerical simulation of the original equations of motion for 3.71 ■ 10^ integration 
time. The parameters are taken from Table 2, 3rd column. The numerically found drift in action 
space is 5.629 ■ 10^"^ and it is by a factor 2 smaller than the upper bound obtained through the 
analytical estimate, which is equal to 1.146 • 10^^. 



stability times in Tables 2 and 3 are of the same order of magnitude, since the values of /i and 
e are comparable at all orders. The stability estimates are checked for the parameters given 
in Table 3 at order 3 (whose stability time is compatible with the computer execution time) 
by comparison with a numerical simulation as shown in Figure 2. The integration time was 
set to be of the order of the stability time (i.e., we integrated up to 3.71 ■ 10^); we found that 
the deviation of the action is bounded as 5.629 ■ 10""^, while the analytical estimate provides 
1.146 ■ 10"'^ (the numerical deviation is therefore bounded with a safety factor 2). 

5.4 A system with oscillating energy 

We conclude by providing an example of a differential system which admits oscillating energy. 
To be more precise, we consider the differential equations 

X = y 

y = —e (sin(x — t) + sin(x)) — sin(x) — rj). 
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The Hamiltonian function for = (in the extended phase space) reads as 

T,x,t) — ^ — e (cos(a; — t) + cos(x)) + T , 

where T is the conjugated action to the time t. For 7^ we get that the variation of the energy 
is given by 

dE dH . dH . dH- dH . o . , , 

^ ^ + + ^ + = '^"^^^ ' ^^^^ • 

Since the normal form equations will provide that = 0, we can conclude that the energy is 
oscillating. The normal form solution to second order provides the following expressions for 
the transformations: 

sin(x — t) sin(5;) 



^2,{y,x,t) 2{y - l)y{2y - I) 8(^-1)3 2{y - l)y 8y^ 

MY,X,t) = cos(X) 



y- 1 y 

sin(2x — t) sin(2x — 2t) sm{t) sin(2a;) 



aoiiY,X,t) 
MY,X,t) 



sin(X) 
Y 

cos(2X) 
2Y 



ycos(2X-t) (F-2)cos(t) cos(2X) 
sin(2X) 



ao2{Y,X,t) 



4y2 



^ sin(2X-t) (F-3)sin(t) sin(2X) 

' ' ^ 2(y-i)(2r-i)2 ^ 2(r-i)2 sy^ ' 

The corresponding normal form at second order is given by: 

^ - ^+ 2(y-i)3y ^' + w-^+o^ie,,) 

Y = 02(6,1,). 

We compare the solution of the normal form equations with the numerical solution by com- 
puting the error as in (42) for £0 = 10"^, /^.q = 10"^ and Yq = ^ (\/5 + l), Xq = 0. Figure 3 
(left panel) shows the error for the 1st (continuous line), 3rd (dashed line) and 5th (dotted line) 
order respectively, up to time IO^tt; the oscillating behavior of the energy is given in the right 
panel of Figure 3. The numerical solution was again obtained using a 4th order Runge-Kutta 
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Figure 3: Left: relative error between the normal form and the numerical solution. We set 
e = 10^^, /U = 10"'^ and 10 = 1 ("\/5 + l); the evolution is in good agreement with the 
solution obtained from the normal form equations (continuous 1st, dashed 3rd, dotted 5th 
order normal form). Right: behavior of the energy of the system, which oscillates around a 
mean value with period 3.86 . 



integration scheme with fixed step size 5t = 10^^. We remark that the difference between 
the numerical and analytical solutions decreases as the order of the normal form increases. 
We conclude by mentioning that the bounds on the small parameters as well as the stability 
estimates can be determined as in Section 5.3. 



6 Appendix A 

We discuss the conditions which must be satisfied by the parameters e, fi, so that the trans- 
formation from original to intermediate variables can be inverted, as well as that from inter- 
mediate to final variables; moreover, we give conditions on the parameters so that the non- 
resonance conditions in the intermediate and final variables are satisfied. Such results rely on 
the following two lemmas which are proven in [29] . 

Lemma A.l. Let yo G M^, {x,t) G T^"''-'^, ro, sq, 5o (6o < sq) be strictly positive parameters 
and let g be a vector function holomorphic on the domain D{yQ, tq, Sq) = {{y, x, t) G C^^+^ : 
\\y — UoW < ■'^0 ; niaxi<j<^ l^l^^j)! < > 1^(^)1 < Sq}. Let us consider the equation 

x' = X + g{y,x,t) ; (43) 

if 

C\\9\L,soe'"'S^' < 1 (44) 
for some positive constant C, then (43) can be inverted as 

X = x + G{y, x', t) , 
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for a suitable function G such that 

\\^\\ro,so-So ^ llfl'llrojSo ■ 

Lemma A.2. Let yo G W\ {x, t) e T^+^ and let ro, sq, fo be stnctly positive parameters with 
Tq < ro; let g be a vector function holomorphic on the domain D{yo, ro, sq) = {(y, x, t) e 
£12^+1 . _ ^^11 < ^ maxi<j<^ < ^0 , 1^(^)1 < So}- -Let us consider the equation 

y' = y + 9{y,x,t); (45) 

if 

<^ll^l|r-o,.o^^ < 1 (46) 

To - ro 

for some positive constant C, then (45) can be inverted as 

y^y' + G{y',x,t) , 

for a suitable function G such that 

\\G\\r'^,so < \\9\\ro,so ■ 

We remark that a careful evaluation of the constant C in (44) and (46) shows that it can be 
fixed as C = 70. 

6.1 Inversion of the conservative transformation 

Let us recall the transformation (7) as 

X = x + 'il)f\y,x,t) 

y = y + V'fny,^,^), (47) 

that we wish to invert as 

X = x + r(^'^)(y,x,i) 

y = y + V^y'''\y,x,t). (48) 

Let fo < ro, 5o < Sq, sq = Sq — 5q; the inversion of the first in (47) can be performed provided 
that 
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with 

l|r(^'^^ll.-o,.-o<IIV'rilfo,.o- 

The second in (48) is obtained from 

y = y + ^ljf){y,X + T^^'^) {^y , x,t) ,t) ^ y + T^^^'^) {y, ~X, t) , 

where 

l|ro,so -i \\Yx \\ro,so ^ WYxx lko,soll^ Ikcso ■ 

6.2 Non-resonance condition after the conservative normal form 

Taking into account (3), we want that the non-resonance condition is satisfied in the interme- 
diate variables, say for a > 0: 

\uj{y) ■k + m\>'^ , \k\ + \m\<K, (49) 

where from (47) we get 

y = y + (50) 
for a suitable function R^^^. In fact, the second of (47) can be inverted as in (50) provided 

for Tq < fo with 
Then we have 



\jj{y) • A; + m| > \io{y) ■k^m\ - eK\R^^\fi^^s^\ujy%^ >a- 



a a 
2 ~ 2 ' 



provided 

e < 



2K\\R(^)\\r'^,J\uj, 



y\\ro 



6.3 Inversion of the dissipative transformation 

Let us now discuss the inversion of (8) that we write for short as 



X = x + a^^\y,x,t) 

Y = ^ + (51) 
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we invert (51) as 

y = y + A(^''^)(y,x,t), 

provided e, are sufficiently small. In fact, the first of (51) can be inverted provided 

where < sq. Inverting the equation as 

x = X + , 

we have 

11^ llfo,So-5o — II" lko,^io • 

Writing the second of (51) as 

y = ^ + /3 W X + ^(-'^) {y^X,t),t)^y + B^y'''^ {y, X, t) , 
we can invert it as 

y^Y + A^y'''\Y,X,t), 

provided 

To — rio 

with So < §0 — So, Ro < ro, being 

||A(-'^)||^,5o< II 11.0,50- 

Notice that A^^'^^ can be bounded as 

11^^^'"^^ 11.-0,50 < ||/3(^)||.o,.-o + ll^f ^ll.o,.-oll>l^^'^^ll.-o.5o ■ 

6.4 Non-resonance condition after the dissipative normal form 

We now turn to the fulfillment of the non-resonant condition in the new set of variables 

|a;(y) • A; + m| > , |A;| + |m| < K . 
To this end, we use the transformation 

and using (49) one easily finds 

\uj{Y) • A; + m| > \uj{y) ■ k + m\ - K\\ujy\\ro\\P^^^\\fo,so > ^ ~ ^ ' 
provided the following smallness condition on the parameters is satisfied: 

i^KILII^^^^Ik.o<^- 
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7 Appendix B 

Lemma B.l. Let f = f{y, x, t) be an analytic function on the domain Cro(^) x C'soC^^^^)- 
Letf>''{y,x,t) = j:ij,m)ez^+^,\j\+\m\>Kfjm{y) e^(^'-+-*) and let < ao < Sq- Then, there 
exists a constant Ca = C'a((To, K), such that 

Wr^'lUso < C„||/||.„,„+.,e-(^+i)'^» , (52) 

with 

a.e(--^)^fi±^j . (53) 
Proof. From the properties of analytic functions, one has that 

Therefore one finds 

(j,m)ez'^+i,|i|+|m|>ii: 

< ll/II^Q^Q+o-g ^ g-(so+ao)(|i| + |m|)g3o(|j| + |m|) 

(i,m)GZ';+i,|i|+|m|>X 

= ll/ll.o,.o4-.o E e--('^l+H). 

(j,m)GZ«+i,|i| + |m|>X 

Taking into account that 

e-M\j\+\m\) < ^-(K+i)^ e-^(bl+H) 

{j,m)eZ(+\\j\+\m\>K {j,m)eZ(+\\j\+\m\>K 



1 — e 2 



one obtains (52) with Ca as in (53). □ 
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